Normal persistent currents in proximity-effect bilayers 



O 

(N 

< 



c 

o 
o 

Oh! 

•4-* 

o3 



o 
o 



> 

On 
O 



X 



0. Entin-Wholman, 1 ' * H. Bary-Soroker, 2 A. Aharony, 1, t Y. Imry, 3 and J. G. E. Harris 4,5 

1 Department of Physics and the Use Katz Center for Meso- and Nano-Scale Science and Technology, 
Ben Gurion University, Beer Sheva 84105, Israel 
Department of Physics, Ben Gurion University, Beer Sheva 84-105, Israel 
:i Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel 
^Department of Physics, Yale University, New Haven, CT 06520, USA 
J Department of Applied Physics, Yale University, New Haven, CT 06520, USA 

(Dated: August 12, 2011) 

We calculate the contribution of superconducting fluctuations to the mesoscopic persistent current 
of an ensemble of rings, each made of a superconducting layer in contact with a normal one, in the 
Cooper limit. The superconducting transition temperature of the bilayer decays very quickly with 
the increase of the relative width of the normal layer. In contrast, when the Thouless energy is larger 
than the temperature then the suppression of the persistent current with the increase of this relative 
width is much slower than that of the transition temperature. This effect is similar to that predicted 
for magnetic impurities, although the proximity effect considered here results in pair-weakening as 
opposed to pair-breaking. 



I. INTRODUCTION 

The average persistent current 1 ' 2 of a large number of 
mesoscopic metallic rings can be used to deduce the sign 
and the magnitude of electron-electron interactions in the 
metal forming the rings. The size of the average current 
is expected to increase with the strength of the interac- 
tions, and its sign reflects the nature of the interactions: 
the magnetic response at low flux is paramagnetic (dia- 
magnetic) when the electronic interactions are repulsive 
(attractive). 3,4 For a large ensemble of rings, the current 
is expected to be periodic in the magnetic flux, with the 
period corresponding to one half of the flux quantum, 
h/2e. The theoretical analysis of Refs. 3 and 4 was 
motivated in part by early measurements of the average 
persistent current in an array of 10 6 copper rings, 5 whose 
sign and magnitude could not be accounted for by nonin- 
teracting electrons alone and therefore should be affected 
by electronic interactions. These experiments confirmed 
the above periodicity, also suggesting that the average 
magnetic response is induced by interactions. 3,4 Similar 
results were later observed on an array of 10 5 GaAs rings 6 
and on an array of 10 5 silver rings. In contrast, measure- 
ments on single rings 8-11 showed the h/e periodicity. In 
an array of 30 gold rings 12 both the h/2e and the h/e 
harmonics were observed. In this paper the authors were 
unable to say whether the h/2e signal was the second har- 
monic of the typical contribution or the first harmonic of 
an average contribution. Overall, the sign of the h/2e 
harmonic measured on metallic rings seems to indicate 
that the low-flux response is diamagnetic, 7,12 implying 
attractive interactions. Recently, Bleszynski-Jayich et 
al. 13 found that the average current in aluminum rings, 
subject to high magnetic fields, is negligible, but typ- 
ical mesoscopic fluctuations remain almost unaffected. 
It seems that these latter experiments can be explained 
within the framework of noninteracting electrons. 14,15 

Interestingly enough, it turned out that the bona fide 
values of the attractive interactions required to explain 



the persistent-current data of the copper 5 ensemble for 
example, would have implied that this metal is supercon- 
ducting at measurable temperatures, of the order of lmK. 
In fact, early experiments on the magnetic response 16 
and on the thermal conductivity 17 of proximity-effect 
systems, whose normal parts were copper and silver, 
also indicated a minute attractive interaction in these 
metals. 18 However, these early measurements allowed 
a broad range for the magnitude of this interaction, 
and therefore did not open a discussion of the reasons 
for the absence of superconductivity in experiments on 
these metals. The latter puzzle became obvious only af- 
ter the measurements of the persistent current on cop- 
per, which requires a transition temperature of lmk. 
Superconductivity has not been detected also in gold 
and silver, and this fact has remained unexplained for 
many years. A possible explanation for this apparent 
puzzle 19 was offered in Refs. 20 and 21, which ar- 
gued (for the first time) that the existence of (seem- 
ingly unavoidable 22 ) tiny amounts of magnetic impuri- 
ties may detrimentally affect superconductivity in such 
metals, reducing their transition temperatures to unde- 
tectable, even zero, values, while leaving the persistent 
current almost unharmed. This stems from the disparity 
of the energy scales determining the renormalized elec- 
tronic interaction pertaining to each phenomenon. The 
interaction-induced persistent current is proportional to 
the renormalized interaction on the scale of the Thouless 
energy, E c = HD /L 2 (where D is the diffusion coefficient 
and L is the circumference of the ring) . Superconductiv- 
ity is lost, however, when the spin-flip rate of the mag- 
netic impurities, h/r s (in units of energy), becomes com- 
parable to the bare transition temperature of the material 
(in the absence of any pair-breaking or pair-weakening 
agents), T^ . In other words, the actual superconducting 
transition temperature T c is determined by the renor- 
malized interaction, on the scale of max[T^, h/r 3 ]. It 
follows that a concentration of magnetic impurities such 
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that 

k B 1% < K/t s < E c , (1) 

will hardly affect the magnitude of the persistent current, 
concomitantly suppressing the superconducting transi- 
tion temperature (below we often use units in which 
H = ks — 1)- Indeed, detailed analysis 21 of the persistent 
current data reported in Refs. 5 and 7 led to the conclu- 
sion that T c g of copper (gold) is in the mK (a fraction of 
mK) range. 
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FIG. 1: (color online) Proximity-effect ring in which the width 
of the strong (weak) superconductor is ds (rfjv), threaded by 
a magnetic flux $ (measured in units of the flux quantum). 

This theoretical picture can be tested, for instance, by 
investigating rings made of known low-superconducting- 
transition-temperature materials, in which a controlled 
concentration of pair breakers have been added. 13 It has 
also been noted that the magnetic flux itself acts as a pair 
breaker, causing a periodic decrease of the transition tem- 
perature but a lesser decrease in the persistent current. 23 
It is interesting to check whether there exist other situa- 
tions where the superconducting transition temperature 
is lowered by some pair-breaking or pair- weakening mech- 
anism, but the (superconducting fluctuation-induced) 
persistent current remains large far above this transi- 
tion temperature. In the present paper we consider this 
question for superconducting-normal (SN) bilayers, e.g. 
made of Al and Cu. 19 Bilayers made of Al and Ag might 
even be better, as they avoid magnetic impurities. The 
'normal' metal could also be a weaker superconductor, 
with a lower transition temperature. The proximity ef- 
fect is known to cause a decrease of the transition tem- 
perature of the bilayer with the relative thickness of the 
S layer, 24,25 and it is interesting to find out what hap- 
pens to the persistent current, which is induced by super- 
conducting fluctuations. This possibility is in particular 
intriguing: unlike the magnetic impurities, the proximity 
effect is not a bona fide pair-breaker, since time-reversal 
invariance is not broken by it. The proximity effect just 
leads to pair-weakening, by 'diluting' the superconduct- 
ing fraction. 26 

Here we present a calculation of the disorder-averaged 
persistent current averaged over an ensemble of bilayer 



rings, each having the geometry depicted in Fig. 1. These 
rings consist of two adjacent metallic rings, with differ- 
ent transition temperatures. The area inside the rings is 
penetrated by a magnetic flux $, which is measured in 
units of the flux quantum hc/e. Below we use the sub- 
script S for quantities characterizing the layer with the 
higher transition temperature, and the subscript TV for 
the quantities belonging to the other one, which may or 
may not be a superconductor. For simplicity, we con- 
fine our calculation to bilayers in the Cooper limit: 24 
this limit is reached when the width of each of the lay- 
ers, ds or djsi, is much smaller than the respective co- 
herence length. 27 Our aim is to explore the possibility 
to deduce the scale of the renormalized electronic in- 
teraction by analyzing concomitantly the superconduct- 
ing transition temperature and the (superconducting) 
fluctuation-induced average persistent current. In other 
words, we examine the persistent current as a function 
of the N— slab relative thickness, and find parameter 
regimes where it is affected much less than the transi- 
tion temperature. 

Since pair breakers, notably magnetic impurities, seem 
to be ubiquitous in several of the metals used in the 
persistent-current measurements, it is interesting to in- 
vestigate their effect in a proximity-effect configuration. 
For instance, it is plausible that in Al/Cu rings, the cop- 
per (the N— slab in our notations) may well include a 
tiny amount of magnetic impurities. We therefore in- 
clude scattering off such impurities in our expressions. 

The transition temperature of an SN proximity bilayer 
in the the Cooper limit is known 24 ' 25 to be determined 
by the effective (dimensionless) electronic coupling, X NS , 
which is the weighted sum of the effective couplings of the 
separate slabs, As (which is positive, since the S— slab is 
superconducting) and Ajv (which may take both signs): 

^ns = Pn^n + Ps^S i (2) 

with 

PN{S) = ^V(S)-^jv(S)M/cff > (3) 

where N"n(S) denotes the density of states at the Fermi 
energy per unit length of the normal (superconducting) 
layer, and 

K s =Af N d N +Af s d s . (4) 

The mean- field transition temperature, , of the bi- 
layer (without magnetic impurities) is then given by 24 ' 25 

where \& is the digamma function whose asymptotic ex- 
pansion, valid for large arguments, is given by 

lim *(z) -►lni: . (6) 

z— >oo 
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The Debye frequency ui D in Eq. (5) (assumed to be iden- 
tical for both slabs) marks the upper cutoff on the effec- 
tive interactions. The result (5) is obtained assuming 
that the two layers are in a good electrical contact. (The 
effect of a barrier between the two layers has been con- 
sidered by McMillan. 28 ) In the limit lo d > T^ s one may 
use Eq. (6) to obtain 
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(7) 



where T^ is the bulk transition temperature of the clean 
S— slab. When the N— slab is also superconducting 
(i.e. X N > 0), T^ s remains finite for all p N (al- 
though quite small for large p N and small X N ). How- 
ever, when X N < 0, the transition temperature of the 
bilayer T^ s approaches zero at a quantum critical point, 
p N = 1/(1 — X N /X S ). The approach is exponential, with 
zero slope (see Fig. 2). In practice, T^ s becomes very 
small for p N > 1/2. In some sense, this inequality re- 
places the left hand side of Eq. (1). As we show below, 
the persistent current remains rather large even in this 
regime. 

The effect of genuine pair-breaking mechanisms on the 
transition temperature was considered a long time ago. 
The seminal paper of Abrikosov and Gorkov 29 found that 
the transition temperature T^ is reduced by magnetic 
impurities to T c s , 
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(8) 



where s = l/(2itT^T s ). This expression is shown in 
the inset in Fig. 2. Unlike T^ s , T c s approaches zero 
at s = exp[*(l/2)] = l/{4r/ E ) » 0.140365 (namely at 
H/t s ~ 0.9T^ ), with a finite slope. Here, j E is the Eu- 
ler constant. This difference in slope between the two 
mechanisms probably reflects the difference between pair- 
weakening and pair-breaking. 26 

In a complete analogy with Eq. (8), a small amount of 
pair-breaking impurities in the N~ slab lowers the tran- 
sition temperature of the sandwich from T^ s to T^ s , 
given by 25 



1 



X 



A ! A' 



1 



2 2nT" s 



1 



2 2-kT? s t. 



Pn 



(9) 



generalizing Eq. (5). 

The rest of this paper describes the calculation of the 
average persistent current, pertaining to a large ensemble 
of bilayers. Section II outlines the derivation of the ef- 
fective Ginzburg-Landau theory for this case, with some 
technical details given in Appendix A. Some quantitative 
results are presented in Sec. III. Since the fluctuations 
are calculated within the high-temperature Gaussian ap- 
proximation, which is valid only above the Ginzburg crit- 
ical regime, Sec. IV presents a critical discussion of this 
regime. That section also contains our conclusions. 
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FIG. 2: (color online) The transition temperature of a 
proximity-effect sandwich in the Cooper limit [Eq. (7)] as 
a function of the weighed width of the normal layer for 
Xn = 0, 0.05, and —0.05 (solid, dashed, and dotted curve, 
respectively). In all curves As = 0.17. The inset depicts the 
Abrikosov-Gorkov expression, Eq. (8), for the reduction of 
the transition temperature by pair breakers. 



II. THE PERSISTENT CURRENT OF A 
PROXIMITY-EFFECT SANDWICH 

Here we present a microscopic derivation of the free en- 
ergy which determines the superconducting fluctuations. 
The Hamiltonian of the bilayer is similar to that used in 
Refs. 20 and 21, 



dvU{v) , 



(10) 



with 



H(r)=^v4WHi CT l(r)Vv(r) 

£T£T' 

-U(r)^4(r)V'J(r)VhX r )V ; t( r ) , 



(11) 



where V'i-C 1 ') creates an electron with spin a at r. The 
interaction V(r) depends on the spatial coordinate x (see 
Fig. 1), 



V(x) = X s /Af s , 
V(x) = X N /Af N , 



-dg <X<0 

< x < d N . 



(12) 



(Note that here the APs are the densities of states per 
unit volume of the two layers.) The single-particle part 
of the Hamiltonian (11) reads 



<(r) 



where 



H Q = [-iV + (e/c)A(r)] 2 /(2m) - n 



(13) 



(14) 



and n is the chemical potential. With the choice A = 
B x r/2, the vector potential A points along the circum- 
ference of the ring in the anticlockwise direction. The 
disorder potential is u(r) = u 1 + u 2 er ■ S, yielding scat- 
tering off nonmagnetic impurities (scaled by %) as well 
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as off magnetic impurities (scaled by u 2 , S denotes the 
magnetic impurity spins. The impurities are modeled by 
point-like scatterers 30 ). 

The quantum partition function Z is 31 



where the action S is 

(16) 

and /3 = 1/T. Here, the annihilation and creation field 
operators in the Hamiltonian (11) (ip and tjj') are re- 
placed by the spinor Grassmann variables ip{r,r) and 
■0(r, r), respectively. 

Since the calculation of the partition function is rather 
technical, we present it in Appendix A. We first 
perform this analysis in the absence of the magnetic 
flux. A Hubbard-Stratonovich transformation replaces 
the Grassmann variables ip and ip by complex bosonic 
variables A(r,r) and A*(r, r) (which are now functions 
of the imaginary time r), and the action is expanded in 
powers of these variables. The Gaussian approximation 
uses only the quadratic terms in this expansion. In the 
Cooper limit, the Fourier transformed bosonic variables 
take only two values as a function of x, namely A s (q, v) 
for — d s < x < and A N (q, v) for < x < d N , where 
q is a two-dimensional vector perpendicular to x. The 
quadratic action then becomes 



q v 

- c[A^(q, i/)A s (q, v) + c.C.]) , (17) 
where Af e s is given in Eq. (4) , and 

P N(S) 2 MCI 

a N( S ) = ~\ Pn(S)1 > c = PsPnI ■ (18) 

A N(S) 

Here, P N r S -> are given in Eqs. (3). The function 7(q, v, T) 
is given by 



7(q 1 i/ ) T) = irr^7(q ) i/,w) , 



(19) 



where 



7(q, v, U ) = (M + M + ^ + \ DcS c?) 1 , (20) 

uj = 7rT(2m + 1) and v = 2ttT£ (with integer m and 
£) are the fermionic and bosonic Matsubara frequencies 
and D eS is the effective diffusion coefficient of the double 
layer, 



Since the sum over oj in Eq. (19) is cut-off by the Debyc 
frequency uj d one finds 



Z = j P«(r,T), Hv, t)} cxp[-S] , (15) +(p / 



1 |z/|+D cff q 2 + 2p N /T 8 + 2bJ D 



4ttT 



(22) 



The bilinear form in Eq. (17) is diagonalized by the 
transformation (for convenience, we omit the explicit no- 
tations of q, v, and T in part of the expressions below) 

A N = u_A_ +u + A + , A s = w + A_ - u_A + , (23) 

where 



= (1 : "v <<<v 



4k 



(24) 



One then finds 

s 2 =w cS EE( a - ^ T ) i A - ") i 2 

q iy 

+ a + (q^,r)|A + (q,^)| 2 ) , (25) 

with 

a ± = (a N + a s )/2 ± k . (26) 

Within this Ginzburg-Landau-like model, the phase 
transition occurs when the first coefficient a ± (q, v, T) 
vanishes as the temperature T is lowered. Since a + — 
a_ = 2k > 0, this transition happens when a_(0,0,T) = 
(while a + remains positive). Equations (18) and (19) 
imply that 

^\ ns _ 7(q, v, 1 )) . 

(27) 



a + a_ = a N a s — c = 



Therefore, at zero flux the transition occurs at T^ s 
which obeys the equation A^g = 7(0, 0, T^ s ). Using 
Eq. (22), this reproduces Eq. (9). 

Finally, we incorporate the magnetic flux into the ex- 
pressions for the action and for the partition function. 
To lowest order (neglecting the effect of the field on the 
order parameter) it suffices to replace q by 32 



q -> q + (2e/c)A . 



(28) 



This follows directly from Eq. (14), remembering that 
the momentum q relates to a bosonic Cooper pair. For 
the circular geometry at hand, the component of q along 
the ring circumference, q„ , becomes 



D eS = PN D N + PS D S ■ 



(21) 



27T 

g y = _(„ + 2*) 



(29) 
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with integer n. The transition is then shifted, with 23 
1 



A 



NS 



= 7 (<7|| =47r$/L,0,T c Ars ($)) 



(30) 



and the persistent current is given by 
^STlnZ 

Within this Gaussian approximation, the fluctuations 
contribution to the partition function can be obtained 
straightforwardly. One finds 



nn 

q v 

nn 



i 



o+(q, u,T)a_(q,u,T) 
1 



N S 



(32) 



where (flux- and temperature-independent) multiplica- 
tive factors have been omitted. Interestingly, this ex- 
pression for the partition function has exactly the same 
form as that found for the 'superconducting' ring in Ref. 
21. The only modification is that now A^ 1 is replaced 
by A^yg. The following calculations thus use the same 
calculational techniques employed in that reference. 



where 



x , , = irriA 

polc/zcro \ 



T 



2tt£L 



1 + 1*1 



Pn 
irTTc 



2F 3 , , 

polc/zcro 



1/2 



with F£„,„ = —j and with F| ( 



(37) 

being the solution of 



*(FiJ)=ln[T™/(4 lE T)}. 

We now present several plots of the persistent current 
based on Eq. (36). Figure 3 shows the first harmonic 
of the current (divided by its value at p N — 0) as a 
function of pjy, for = 1/t s = 0. To avoid critical 
fluctuations (see below), we restrict ourselves to a rela- 
tively high temperature, T = 4T c q. The same figure also 
shows the transition temperature for the bilayer, divided 
by T^ . Clearly, the relative persistent current decreases 
much more slowly than the relative transition tempera- 
ture. This slower decrease is similar to that found in Refs. 
20 and 21, resulting from the effects of pair breakers. As 
an example, for the parameters used in Fig. 3, the tran- 
sition temperature at p N = 0.7 is very small, T^ s {p N = 
0.7) ~ 10~ 6 Tj) while the first harmonic of the current is 
given by I{p N = 0.7) = 0.16/ [p N = 0) = 0.13£: c . 



III. RESULTS 

Since the important contribution to the persistent cur- 
rent comes from the zero transverse mode (perpendicular 
to the x— direction), 15 ' 20,21 we replace the sum over q by 
a one- dimensional summation over the discrete values 
of q,,, Eq. (29). Assuming that the Debye frequency uj d 
is the largest energy in the problem, the denominator in 
Eq. (32) becomes 21 



A^-7(q,^,T) 



In 
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F(n,£) 



1 



and therefore the persistent current is 21 

(n + 2$)*'(F) 



ln(T/T» s ) + *(F) - *(l/2) 



(33) 



(34) 



(35) 



(We remind the reader that E c — D cff /L 2 is the Thou- 
less energy). As shown in Ref. 21, this expression for 
the persistent current can also be written as a Poisson 
summation, 



I = — 4eT sin(47TTO$) 

oo 

X [ eX P( 27r * a; zerc 
I j=l 



exp(27rixi j ole )] , (36) 
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FIG. 3: (color online) The transition temperature of the bi- 
layer T^ s , normalized by T^ (blue solid line) and the first 
harmonic of the current, divided by its value at p N — 0, from 
Eq. (36) (red dashed line). The green dash-dotted line shows 
the approximation (41) for the current. The parameters used 
for all graphs are T|, = 1.27K (equivalent to A s = 0.17), 
uj D = 400K, T% = (i.e. X N = 0) and l/r s = 0. The current 
is plotted for T = and E c = lOTj,. 



At a fixed p N , the persistent current decreases with 
increasing temperature. Figure 4 shows the current (in 
units of the Thouless energy E c ) as a function of the 
temperature for a specific choice of the parameters and 
for three values of p N . Each of these plots shows the 
current only above the transition temperature T^ s . As 
anticipated in the Introduction, the persistent current 
increases with increasing E c . This can be seen from Eqs. 
(36) and (37), in which the decay of I is determined by 
the ratio T/E c . Using the relation j\f cS T^ = g/[E c /T^], 
where g is the dimensionless conductance, the parameters 
used in Fig. 4 are equivalent to g = 1000. 

Finally, we discuss the effects of a positive transition 
temperature and a finite amount of magnetic impu- 



6 




FIG. 5: (color online) Same as Fig. 4, but with T% = O.OOlif 
(i.e. X N — 0.077). The current is plotted for p N — 0.7 (blue 
dashed) and p N = 0.8 (red dash-dotted). Thick (thin) lines 
represent 1/t s = (l/r s = E c ). The graphs are plotted for 
T > T* s . 



FIG. 4: (color online) The first harmonic of the current, in 
units of the Thouless energy E c , versus T/E c , for T^ = 
1.27K, E c = 0.015K, lo d = 400K, T% = l/r s = and 



Af c ffTc = 10 J . The current is plotted for p N = 0.45, 0.5, 
and 0.7 and for T > T$ s (p N ). 



rities in the normal slab. Figure 5 shows the first har- 
monic of the persistent current versus the temperature 
for = lmK, which is the estimated minimal value for 
the pure transition temperature of copper derived in Ref. 
21, with and without magnetic impurities. We see that 
at high temperatures the current is not very sensitive to 
the pair breaking. As might be expected, the weak su- 
perconductivity of the N layer causes an increase in the 
persistent current. 

Interestingly, both Fig. 4 and Fig. 5 exhibit 
fluctuation-induced persistent currents which are much 
larger than the Thouless energy E c , at temperatures 
above the superconducting transition temperature of the 
S material. This persistent current increases, and its de- 
cay with temperature becomes slower, as E c increases. 
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The above plots were based on Eq. (36), which sums 
over many values of i and j. For example, Figs. 4 and 
5 used j,£ < 10 4 . We next describe approximate ex- 
pressions, which are valid at high temperatures, when 
T 3> T^ s . We restrict this discussion to the case 
l/r s = 0. Defining the small parameter 



w 



and noting that i£(w- 
and therefore 



-l/ln[T c ^/(4 7£ T)] 

a — 1 /w, we find F£ f 



(38) 

z, m ' J -j+W, 



C 

pole 



T 



2nE c [l + 2 3 + \£\] 



(39) 



Inserting this approximation into Eq. (36), the Poisson 
summation form of the current becomes 



I « -4eT sin(47rm$) ^ exp[-m v / 27rT[i + 2j + \£\}/E c exp(mwy/2wT/Ej\l + 2j + \£\}) - 1 



l 3=0 



—AeTw sin(47TTO$) 

J2 m V^T/E c /[l + 2j + \£\] exp[-mv/27rT(l + 2j + \£\)/E c ] , 

m=l £ j=Q 



(40) 



r 



where the last approximation applies only for wiw\/2'kT jE c <C 1. At intermediate temperatures, 
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when both this condition and w <C 1 are obeyed, the 
current decays as I ~ wl 1 , where I 1 is independent of 
T^ s . Substituting Eq. (7) for T^ s in the expression for 
w, we finally end up with 



ln[4 7£ T/T c s ] 



h HA lE T/T^] + 



Pn( 1 - x n/ x s) 
Pn x n+Ps x s 



(41) 



where I is the persistent current (for the same flux) at 
p N = 0. As seen in Fig. 3, this approximation is quite 
good. 

Unlike the case of the magnetic impurities, in which the 
persistent current remains non-zero even when the tran- 
sition temperature vanishes, in the case of the bilayer 
the persistent current vanishes when T^ s = 0. When 
X N = 0, the transition temperature approaches zero as 
p N increases towards the quantum critical point, which 
occurs at p N = 1. When X N < 0, this critical point oc- 
curs at a threshold p^ < 1, and the fluctuation- induced 
persistent current vanishes above this threshold. How- 
ever, as p N approaches this critical threshold, the current 
decreases linearly with p N [as seen from Eq. (41)]. Since 
the transition temperature decays exponentially towards 
that point, we again find that the persistent current re- 
mains significant even when the transition temperature 
is negligibly small! 



IV. DISCUSSION 

The calculations above were carried out within the 
Gaussian approximation. This approximation usually 
breaks down close to the phase transition, where higher 
powers of the order parameters must be taken into ac- 
count in the expansion of the action S. This happens be- 
low the so-called Ginzburg temperature, T Gi . Therefore, 
one should not trust the above results for temperatures 
in the range < T < T Gi . In this range, we need to 
supplement Eq. (17) by the quartic terms, which should 
be derived by continuing the expansion of Eq. (A2) in 
powers of the A's. In principle, this expansion has the 
form 



x A Q (q 1 ,i/ 1 )A j8 (q 2 ,i/ 2 )A 7 (q3,i/3)A4(q 4 ,i/ 4 ) , (42) 



where a, f3, 7, 8 take the values S or TV, and the sums are 
restricted by £ • Q 4 = Ei v i = °- 

The calculation of the coefficients B a/3rfS goes beyond 
the scope of the present paper. 31 Usually, these coeffi- 
cients are assumed to be independent of the momenta 
q 4 and frequencies za, since such dependencies are less 
relevant near the phase transition in the renormalization 
group sense. Furthermore, Eq. (18) shows that the in- 
teraction terms in Eq. (17), |A S | 2 /A S and \A N \ 2 /X N , 
are equal to their bulk values multiplied by p s and p N , 
respectively. This indicates a renormalization of A s and 



Ajy by the factors p X J 2 and by pj/ 2 ' respectively. In 
analogy, we conjecture that the various coefficients in 
Eq. (42) are also given by their bulk values, multiplied 
by the same renormalization factors. We next replace 
these order parameters by A ± , from Eq. (23). For sim- 
plicity, we restrict the following discussion to the special 
case A 



N 



0. This should suffice to demonstrate our 
arguments. In this special case one has A^-g — Ps^s? 
a_ = a s = p|[A^g — 7], u + = 1 and u_ = 0. Also, 
a + = a N = 00, and therefore we can ignore all the fluc- 
tuations associated with A + . Finally, the quartic action 
becomes 



S 4 = i/W cff B_ ]T ]T f[ A_ (q,,^) 



(43) 



where again J2i — J2i v i = an d we se t B_ = p%B , 
with B Q = 7C(3)/(8tt 2 T 2 ) having the bulk value of the 
quartic term. 33 Keeping only the first term in Eq. (25), 
we find the usual structure of the effective Ginzburg- 
Landau action, except for the renormalization of the co- 
efficients. As we discuss below, the dependence of a_ on 
q and on v is very different from the simple form used in 
standard Ginzburg-Landau theories. 

The literature contains several ways to estimate the 
Ginzburg region. Since here we calculate the persistent 
current, we define that region as the range where the 
Gaussian calculation presented in the previous section 
must be modified by inclusion of the quartic terms. Ex- 
panding the partition function Z to leading order in B_ , 
the free energy becomes 

F = -TlnZ fl + 3AA cffJ B_[^^(|A_(q, i ,)| 2 )] 2 , (44) 



where (. . .) denotes averaging with the Gaussian action, 



(|A_( q ,^)| 2 ) = l/[/3AA cff a_(q^,T)] 



(45) 



The correction to the persistent current due to the quar- 
tic term thus becomes 

J eft v q v q 

The above Gaussian results can be used only if this addi- 
tional contribution is smaller than that calculated above, 
Eq. (35). 

At |$| =0, the denominators in the sums in Eqs. (35) 
and (46) vanish for n = £ = 0, at the critical temperature 
T Jo S ( $ ) which satisfies Eq. (30). Moving slightly away 
from this temperature, i.e. at small T — T^ s (&), each of 
these sums is dominated by its first term, with n = I = 0. 
Most of the discussions in the literature proceed by con- 
sidering only these 'zero-dimensional classical' terms. 33 
Following this 'tradition', i.e. keeping only these leading 
terms in all three sums, and comparing I with 81, we find 
that for $ w the latter can be neglected if 
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(47) 



r 



where s was defined after Eq. (8). The left hand side 
is the denominator for n = £ = 0, which vanishes at 
the mean field transition temperature. Apart from mul- 
tiplicative factors of order unity, one obtains a similar 
'zero-dimensional classical' condition using other defini- 
tions of the Ginzburg region. 33 ' 34 Close to the transition 
at T^ s one usually replaces T by T^ s in the denomi- 
nator of the right hand side. Equation (47) then agrees 
with the usual Ginzburg criterion in d dimensions, which 
would give (T Gl - T c ) ~ \u{T Gl /T c ) ~ T~ 2/(4 ~ d) , with 
d = 0, 33 except for the additional factor l/p% How- 
ever, this substitution is problematic when is very 
small, as in our case. Therefore we prefer to keep T 
also on the right hand side, and then solve Eq. (47) 
as an equality. For p N < 0.5 and for the parameters 
used above, the resulting T Gi turns out to be quite close 
to Ttf) , and therefore the 'zero-dimensional classical' 
Ginzburg region is very narrow. Therefore, Fig. 6 shows 
these two temperatures only for p N > 0.65. As seen 
in this figure, the Ginzburg temperature does not de- 
cay as fast as the transition temperature at large values 
of p N . For example, for the value of Af e sT^ used in 
the figure, we have T Gl (p N = 0.8) = 9 * lO -7 !^, while 
t S) S (pn = 0.8) = 6*10 -11 T^. In any case, the 'classical' 
Ginzburg region is quite narrow. Had we stopped here 
(as done in much of the literature), we would conclude 
that our Gaussian results for the persistent current can 
be used for practically all temperatures (above T^ s ) and 
relative widths of the bilaycr. 



x 10" 




T /T S 

Gi cO 

T NS /T S 

cO CO 



0.65 



0.75 



0.85 



'N 



0.95 1 



FIG. 6: (color online) T Gi /T^ and T$, s /T c s versus p N , for 
Afc«T c s = 10 5 . Note the factor of 10~ 5 on the j/-axis. 



the Gaussian approximation, on which most of our re- 
sults are based, requires the summation over many clas- 
sical and quantum fluctuations. For some reason, most of 
the literature ignores the 'quantum' fluctuations coming 
from non-zero £'s, and keeps only the 'one-dimensional' 
terms with n ^ o. 33 ' 35 As discussed above, and in Ref. 
21, the persistent current is affected by both types of 
fluctuations. The sum in Eq. (35) always converges, 
becoming of order l/ln[T/T^] for T^ «T« E c . Sim- 
ilarly, the last sum in Eq. (46) also converges, becom- 
ing of order l/(\n[T/T^ ]) 2 . In contrast, the first sum 
in Eq. (46) does not converge, and thus it depends on 
the cutoffs imposed on the wave vectors n and on the 
frequencies £. This problem arises since our calculation 
necessitates the replacement of the 'usual' Green func- 
tion 1/(T -T c + v + Dq 2 ) by 1/[A _1 — 7(q, v, T)], with 
a logarithmic dependence at large T, v and q [see e.g. 
Eq. (33)]. Since this sum depends on the cutoffs, the 
resulting Ginzburg criterion will also depend on these 
cutoffs. 36 In our case, the dirty diffusive limit imposes 
the cutoffs \i/\, -D off q 2 <C 1/ti, where r + is the elas- 
tic mean free time. 21 Replacing the sum X)i/q(V a -) by 
some cutoff dependent constant still shows that \SI/I\ 
decreases with increasing T. The details of this cutoff 
dependent criterion go beyond the scope of the present 
paper. 

Our calculation was done for an ensemble of proximity- 
rings in the same plane. Qualitatively, we expect simi- 
lar behavior for two rings which are deposited on top of 
each other, which may be easier to realize experimentally. 
However, the explicit calculation for the latter case still 
needs to be carried out. It is also interesting to calculate 
the persistent current when the bilayer is connected to 
leads. This also remains for future calculations. 

In conclusion, we have demonstrated that the effect 
of pair weakening due to the proximity between a su- 
perconducting ring and a normal (or a weakly supercon- 
ducting) ring is similar to, but not identical with, that of 
pair breaking: the persistent current decays slowly with 
the relative width of the normal layer, and persists even 
when the superconducting transition temperature (which 
decays faster) is very small. Since this relative width can 
be controlled, it would be interesting to check our quan- 
titative predictions experimentally. 



As one moves away from the critical region, each of 
the sums in Eqs. (35) and (46) must be supplemented 
with other terms, involving both non-zero wave vectors 
(n ^ 0) and non-zero Matsubara frequencies (£ ^ 0). 
In fact, we find that as n increases we need to include 
more values of £ to obtain convergence, and that even 
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Appendix A: The partition function 

Applying the Hubbard-Stratonovich transformation to 
Eq. (15), and integrating the fermionic part of the action, 
the partition function is cast into the form 31 



partition function of noninteracting electrons; the second 
one represents the contribution of the superconducting 
fluctuations to that function. Its calculation requires the 
correlation 

/C(r, r', t-t')^ -<Tr{G p (r, r', r - r')a y 

xG h (r',r,r'-T)a y }) , (A7) 

where (...) indicates averaging over the impurity config- 
urations (see Ref. 30 for details). Upon averaging, the 
spatial dependence of JC becomes a function of x, x', and 
p — p', where p _L x, see Fig. 1. Hence, 



Z = J 2?{A(r,r),A*(r,T)} e - 5 , 



(Al) 



TrOnO^- 1 )}! 
dxdx' 



2 „d 



/ ^EE A x'(q^)^(q^)A*(q,^) , (A8) 



with the action 



where 



S = 



dr 



|A(r,r) 

V(x) 



iTrjln^- 1 )} • (A2) K xx ,{q,v) = £ £(Tr{G(x, x', Pl + q, p 2 + q, u + v) 



Pi,p 2 w 



Here f3 is the inverse temperature and Q 1 is the inverse 
Green function at equal positions and imaginary times, 



x (j y G\x',x, -Pi,-p 2 ,-w)a-j,|) , 



(A9) 



Q~ x = 



G~ 



ia y A 



(-ia v )A* G- h 



with 



G~ 



-d T - K 

-1/(0) 



(0) 

tt 



-n 



(0) 

n 



-d r -U 



(0) 

u 



being the particle inverse Green function, and 



Gl 



(o) 
n 



(0) 

tt 



0/(0) 

-d T + H 



(o) 
U J 



(A3) 



(A4) 



(A5) 



and both Green functions are the particle one, 21 i.e., G = 
G p . [In Eq. (A8), d = d N + d s is the total width of the 
sandwich.] We use the notations uj = Lu n = 7rT(2n+l) for 
the fermionic Matsubara frequencies, and v = v t = ttT2£ 
for the bosonic frequencies. Note that q, p x and p 2 are 
two-dimensional vectors normal to x. 

Inserting these results into the expression for the ac- 
tion [sec Eq. (A2)], the Gaussian fluctuation- induced 
partition function, Zr^, takes the form 

Z flj2 = J V{A x ( q ,v),A* x ,( q ,v)}e- s * , (A10) 



with 



being the inverse Green function of the holes. The factor 
/3 was introduced into the last term in Eq. (A2) to keep 
the argument of the log dimensionless (it does not affect 
any of the following discussion) . 

The integration over the bosonic fields in Eq. (Al) 
is carried out using a stationary-phase analysis 31 of the 
action S. At temperatures above the transition temper- 
ature, this amounts to expanding the second term on the 
right-hand side of Eq. (A2) to second order in A (the 
first-order contribution to the expansion being zero) 



5 2 = EE/ dxdx'A* x {^v)[ 



8(x - x') 
V(x) 



-plC xx >((i, v))A x ,(q, v) 



(All) 



(here V is the attractive interaction, in units of energy 
x length). In the Cooper limit, the bosonic variable A 
takes only two values as a function of x, AAr(q, v) for 
< x < djy, and Ag(q, v) for — ds < x < 0. Therefore 
the action becomes 



Tr{ln(^ -1 )} = Tr{ln/3 
rfrdr' f 13 drdr' 



G p 1 





Gl 1 



S 2 = £I>t(q,i/)S(q,i/)A(q,i/) 



(A12) 



n 2 



i3 2 



/C(r,r',T-T')A(r',r')A*(r,T) , 

(A6) 



where f2 denotes the volume of the system (we added the 
factors of volume and (3 to keep S dimensionless). The 
first term on the right-hand side of Eq. (A6) will give the 



where 

S(q,v)=[3 

d 2 N K NN {<\,v) d s d N K NS {ci,v) 







_0 

dgNsXS 1 



1 

d 2 " 



d N d s K SN {q,v) d s K, ss (q,v) 



(A13) 
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7Vjv(s) denotes the density of states of the normal (su- 
perconducting) layer per unit length, and A^(q, u) = 

{A^(q,i/),A£(q,i/)}. 

The functions K. xx >(q,u), Eq. (A9), are calculated by 
extending the method employed in Rcfs. 24 and 25 to 
include the dependence on v and on the two-dimensional 
wave vector q. The calculation is valid in the dirty limit 
[in which (D/2ttT) 1 ^ 2 is much larger than the mean- free 
path of the relevant metal, where D is the diffusion coef- 
ficient] . For simplicity, we omit the vector potential from 
this calculation; its effect is incorporated into the result 
at the end of Sec. II. 

We follow the derivation given in Rcf. 24, and begin 
by presenting the response function K in the form 

1/C xx ,(q,^)=£ff xx ,(q,^). (A14) 

Had the normal part of the bilayer filled the entire space, 
then H xx , = H^} , where 



\2u}\ + \u\ + 2/r s + D N (q L 2 + q 2 ) ' 



We have allowed for scattering off magnetic impurities in 
this metal, whose effect is presented by the spin-flip rate 
1 /t s . (The effect of scattering off nonmagnetic impurities 
is contained in the diffusion coefficient.) As seen from Eq. 
(A15), the function H xx ) obeys a diffusion equation 

(|2w| + \v\ + 2/t s + D N q 2 

-D N ^) H £)( q ,v,u J ) = 2nti N 5(x- X ') . (A16) 

Quite similarly, when the S metal fills the entire space 
one finds 

(|2 W | + M+£ s q 2 

-D s ^)h { x S J{^v^)=2-kM s 5{x-x') . (A17) 

Here it was assumed that the S metal is not doped with 
magnetic impurities. It follows that in order to find H xx > 
of the double layer, one has to solve the set of equations 



(|2w| + M + 2/t s + D N c? - D N -^^H xx ,(q, v, u) = 2nM N 6(x - x') , x' > , 

{\2u:\ + \ V \+D s ci 2 -D s ^)H xx ,(ci, Vi u:) = 2ttN s 8(x~ X ') , x'<0, (A18) 



I 

with the appropriate boundary conditions. Such a where 7 (q, v, w) was defined in Eq. (20). 
scheme has been undertaken in Refs. 24 and 25, lead- 
ing to the result 

II ( \ n ^N -( \ 

H NN (q,v,Lj) = —~ — ~r^(q,v,Lj) , 

d N J\l N + d s N s 

H NS {^v,u) = H SN (q,v,uj) = ^ Tr7 (q,^,cj) , 

d N M N + d s N s 

H ss (q, u, u) = ~ 7 (q, v, to) , (A19) Inserting Eqs. (A14) and (A19) into Eq. (A13) brings 

d N N N + d s Af s the action S, Eq. (A12), into the form (17). 
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